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g ■ Abstract 

■ An asymptotic interation method for solving second-order homogeneous lin- 
^ ■ ear differential equations of the form y" = Xo(x)y' + so{x)y is introduced, where 

■ Xo{x) ^ and so{x) are Coo functions. Applications to Schrodinger type prob- 
lems, including some with highly singular potentials, are presented. 
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1. Introduction 

Second-order homogeneous linear differential equations arise naturally in many 
fields in mathematical physics. There are many techniques available in the liter- 
ature that can be used to solve these types of differential equation with boundary 
conditions. The main task of the present work is to introduce a new technique, which 
we call the Asymptotic Iteration Method, to solve second-order homogeneous linear 
differential equations of the form 



where Xo{x) and so{x) are defined in the some interval, not necessarily bounded, 
and Xo{x) and So{x) have sufiiciently many continuous derivatives. 

2. The Asymptotic Iteration Method 

Consider the homogenous linear second-order differential equation 



where Ao(a;) and so{x) are functions in Co© (o, In order to find a general solution 
to this equation we rely on the symmetric structure of the right hand side of (2.1). 
Indeed, if we differentiate (2.1) with respect to a; , we find that 



y" = Xo(x)y' + soix)y 



(1.1) 



y" = Xo{x)y' + so{x)y 



(2.1) 



y'" = Xi{x)y' + si{x)y 



(2.2) 




y 



,//// 



X2{x)y' + S2{x)y 



(2.3) 




(2.4) 



and 



y 



^^+^^ = Xnix)y' + Sn{x)y 



(2.5) 



respectively, where 




Sn-l + AoA„_i, 



and Sn = s^_i + soAn-i- 



(2.6) 
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From the ratio of the (n + 2)*'* and (n + 1)*'* derivatives, we have 
-^ln(y("+i)) 



dx 



(2.7) 



We now introduce the 'asymptotic' aspect of the method. If we have, for sufficiently- 
large n , 

(2.8) 



Sn-l 



a, 



then (2.7) reduces to 

which yields 

y(^+i)(a;) = Ciexp 



d 



dx Xn-l 



(2.9) 



Xnjt) 

K-i{ty 



( ^ 



dt ] = CiXn-1 exp 



j{a + Xo)dt\ (2.10) 



where Ci is the integration constant, and the right-hand equation follows from 
(2.6) and the definition of a. Substituting (2.10) in (2.4) we obtain the first-order 
differential equation 



y (q! + Xo)dt 



(2.11) 



which, in turn, yields the general solution to (1.1) as 

y{x) = exp ^- J adt^ C2 + C1 J exp ^ J (Ao(t) + 2a{T))dT^ dt 

Consequently, we have proved the following theorem 
Theorem: 

Given Xq and sq in CoQ{a,b), then the differential equation 

y" = Xo{x)y' + so{x)y 
has a general solution (2.12) if for some n > 



(2.12) 
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where 

Afc = A'fe_i + Sk-i + AoAfc-i and Sk = s'^_^ + soAfc-i (2.14) 
for k = 1,2, . . . ,n . 
3. Some illustrative examples 

3.1 Differential equations with constant coefficients 

If so{x) and Xo{x) are constant functions, for example, Aq = 4 and sq = —3 , 
then, for the differential equation y" = Ay' — 3y , the computation of A^ and Sn 
by means of Eq.(2.14) implies 

An = ^(3"+'-l) and = -^(3"+^ - 1). 



The condition (2.13) implies lim ^ = —1 which yields from (2.12) the general 

n— >oo " 

solution y{x) — C^e^ ^ Cie^^, as expected by the application of elementary meth- 
ods. Generally speaking, if we consider the differential equation (1.1) with Ao(a;) 
and So (a;) are constants, we have from (2.14) that 

An — Sn-i + AqA^-i and — sqA^-i. 
Consequently, the ratio -f^ becomes 

An Sn-1 + AoAn-1 + Aq 

^n — 1 

which yields, by means of (2.13), that 

Sn So 



An + Ao 

therefore 

+ ^0^-^0 = 0. (3.1) 

J ^n 

This is a quadratic equation that can be used to find the ratio ^ in terms of Aq 
and So • Therefore, the expected solutions for the differential equation (1.1) with 
constant coefficients follow directly by means of (2.12). 

3.2 Hermite's differential equation 



Many differential equations which are important in applications, such as the equa- 
tions of Hermite, Laguerre, and Bessel, can be solve using the method discussed in 



Asymptotic iteration method for eigenvalue problems 



page 5 



Section 2. As an illustration we discuss here the exact solution of Hermite's Equa- 
tion by means of the iteration method; other differential equations can be solved 
similarly. Hermite's differential equation takes the form 

/" = 2xf' - 2kf, -oo < a; < oo. (3.2) 

Here we have Aq = 2x and sq = —2k . Using (2.14), we can easily show that 

n+l 

6 = Xn+iSn-Sn+iXn = 2''+^llik-i), n = 0, 1, 2, . . . . (3.3) 

Therefore, for the condition 5 = to hold, we must have k a non-negative integer, 
usually known as the order of the Hermite equation. Consequently, for each k , the 
ratio yields 

k = 0, ^ = ^ = ... = 0^/o(x) = l 

Aq Ai 

k=l, ^ = ^ = ...= -- ^ f^{x) = X 

Ai A2 X 



and so on. Clearly, the expressions for the exact solutions fk{x) generate the well- 
known Hermite polynomials. We can easily verify that the general form of fk{x) , 
/c = 0, 1, 2, . . . , is given in terms of the confluent hypergeometric functions [2] by 

f2kix) = i-l)'2'^ (l) iFi f-k- ^; xA , (3.4) 



and 



k 



f2k+i{x) = (-1)'=2'= ( 0^ X ^F^ (^-k; | ) , (3.5) 



where the Pochhammer symbol {a)f~ is defined by (a)o = 1 and {a)k = a{a + 
l)(a -|- 2) . . . (a -|- — 1) for = 1, 2, 3, . . . , and may be expressed in terms of the 
Gamma function by {a)k = r(a -|- /c)/r(a), when a is not a negative integer — m , 
and, in these exceptional cases, {—m)k = if A; > m and otherwise {—m)k = 
{-l)'^m\/{m - k)\. 
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3.3 Harmonic oscillator potential in 1-dimension 

Although, the iteration method discussed in section 2 can be apphed to any second- 
order homogeneous hnear differential equations of the form (1.1) with Aq 7^ , we 
shall concentrate in the rest of the article on the eigenvalue problems of Schrodinger 
type. We shall show that Eq.(2.12) with the conditions (2.13) and (2.14) gives 
a complete solution for many important Schrodinger-type problems. Through a 
concrete example we explore the exact solutions of Schrodinger's equation for the 
harmonic oscillator potentials, namely 



where V ^ L2{—oo, 00) . In the limit of large x , the asymptotic solutions of (3.6) 
can be taken as any power of x times a decreasing Gaussian. With this in mind we 
write the 'unnormalized' wavefunctions as 



where the functions f{x) to be found by means of the iteration procedure. Substi- 
tuting (3.7) into (3.6), one obtains 

which, by comparison with (3.2), yields the exact eigenvalues 

En = 2n+1, n = 0,1,2,... 

and the functions fn{x) , n = 0,1,2, ... , are the Hermite polynomials obtained 
above. Therefore, using (3.4) and (3.5), the unnormalized wavefunctions of the 
Schrodinger equation (3.6) are: 




(3.6) 



'i/j{x) = e 2 f(^x) 



(3.7) 



M^) = (-i)"2' 



( 



2 



1 




(3.9) 



for n = 0, 2, 4, . . . and 



= {-ir2 



in 




n 



xe 2 iFi —n: 



3 2^ 
' 2' " j • 



(3.10) 



for n = 1, 3, . . . . The normalization constant of ip{x) can be computed by means 
of 1 1-01 1 = 1, as we shall shortly show. 



Asymptotic iteration method for eigenvalue problems 



page 7 



3.4 Gol'dman and Krivchenkov Potential 

The Gol'dman and Krivchenkov Hamiltonian is the generahzation of the Harmonic- 
oscillator Hamiltonian in 3-dimensions; namely 

■^^r^ + 2^)^ = E^ (3.11) 

where ip e L2(0, cxo) and satisfies the condition ■0(0) = known as the Dirichlet 
boundary condition. The generalization lies in the parameter 7 ranging over [ 0, 00) 
instead of the angular momentum quantum number I — 0,1,2, ... . For large r , the 
exact solutions of (3.11) are asymptotically equivalent to the exact solutions of the 
harmonic-oscillator problem with eigenfunctions vanishing at the origin. Therefore, 
we may assume that the unnormalized wavefunction if; takes the form 

^(r) = r^+ie-'^/(r) (3.12) 

where, again, /(r) is to be determined through the iteration procedure discussed 
in Section 2. Substituting (3.12) into (3.11), we obtain 

-2(r-2±l)^ + (2^ + 3-E)f. (3.13) 



dr^ \ r J dr 

where Ao(r) = 2(r— -^^) and so{r) = 2^ + 3 — E . By means of Eq.(2.14) we may 
compute Xn{r) and Sn{r). That result, combined with the condition (2.13), yields 

Eo = 3 + 27, El = 7 + 27, ^2 = 11 + 27, ... 

respectively, that means 

= 4n + 27 + 3, for n = 0, 1, 2, . . . . (3.14) 

Furthermore, with the use of /(r) = exp(— J j^dr) , Eq.(2.13), after some straight- 
forward computations, yields 

f (r)-Y^( ilfeon-2fc r(n + l)r(27 + 2n + 2)r(7 + n-fe + l) on-2fc 

"^"^ ^ ' r(/c + l)r(n - /c + l)r(n + 7 + l)r(2n + 27 - 2A; + 2) 

(3.15) 

In order to show that (3.15) together with (3.12) yields the exact wavefunctions for 
the Gol'dman and Krivchenkov Hamiltonian, wc may proceed as follows. Using the 
Pochhammer's identity (a)_fc = r(a — /c)/r(a) , we can write (3.15) as 



f _ 2 (7 + n + l)-fc / J_\ 

Mr)-2r 2J. k\{n+lU{2n + 2^ + 2U,[2r) 



2k 
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Since (a)_fe = (— 1)^^/(1 — a)k we have 



;„M = 2V"til=Mz^^2^(-L)' (3.16) 

fc=o ■ ^ ^ 



in which we have used Gauss's duphcation formula 

a\ /a + 1 



(a)2n = 2 



2n 



2 /„ V 2 



The finite sum in (3.16) is the series representation of the hypergeometric function 
2F0 . Therefore 

U{r) = 2V" 2F0 (-n, -n - 7 - 1; -; -1 
= (-1)'^2'^ n! Ll^Hr^) 
= (-ir2-(^7 + iFi(^-n;7 + |r2 

Consequently, the unnormalized wavefunctions take the form 

^|;{r) = {-1)^2^ (^7 + ^) r^+^e'^iFi (^-n;7+ ^;r^^ 

for n = 0,1,2,.... The normalization constant for ?/'(r) can be found using 
11-011 = 1 which leads to the exact wavefunctions of the Gol'dman and Krivchenkov 
potential, namely 



^''•> = t-^>"'/ n!r(7+i) '-'""'""''''(""' 

Some remarks are in order: 

1. The exact odd and even solutions of the harmonic oscillator potential in 1- 
dimension can be recovered from (3.17) by setting 7 = and 7 = — 1 respec- 
tively. 

2. The exact solutions of the harmonic oscillator potential in 3-dimensions can be 
recovered from (3.17) by setting 7 = where I = 0,1,2,... is the angular 
momentum quantum number. 

3. The exact solutions of the harmonic oscillator potential in -dimensions can 
be recovered from (3.17) by setting 7 = Z -|- ^(A'" — 3) where N >2 . 
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4. Singular Potentials 

We discuss in this section the apphcation of the iteration method discussed 
in section 2 to investigate two important classes of singular potentials. The first 
class characterized by the generalized spiked harmonic oscillator potentials which 
has a singularity at the origin, and the second class characterized by the quartic 
anharmonic oscillator potentials where the perturbative term diverges strongly at 
infinity. Different approaches are usually applied to deal with each of these classes. 
The asymptotic iteration method, however, can be used to investigate the eigenval- 
ues for both classes of potentials. 

4.1 Generalized Spiked Harmonic Oscillator Potentials 

Since the interesting work of Harrell [7] on the ground-state energy of the singular 
Hamiltonian 

(P A 

H= — — + a;2 + — , a; e [0,oo), A > 0, a > 0, (4.1) 

known as the spiked harmonic oscillator Hamiltonian, the volume of research in this 
field has grown rapidly. A variety of techniques have been employed in the study of 
this interesting family of quantum Hamiltonians [4]-[14]. We shall investigate here 
the solutions of the spiked harmonic oscillator Hamiltonian (4.1) in arbitrary dimen- 
sions by means of the iteration method. That is to say, we examine the eigenvalues of 
the Hamiltonian known as the generalized spiked harmonic oscillator Hamiltonian 



, .2 , 7( 7+1) 

x'^ x^ 



H = + '^'^^ ' + —, x e [0,oo), A > 0, a > 0, (4.2) 



where 7 = Z + (A^ — 3)/2 for N > 2 . Clearly we may compute the eigenvalues of 
the spiked harmonic oscillator Hamiltonian (4.1) directly by setting N = 3 and 
1 = in (4.2). 

The wavefunctions (3.12) of the Gol'dman and Krivchenkov Hamiltonian sug- 
gests that the exponent in the power of x term of the exact solutions of (4.2) should 
depend, at least, on the parameters 7 and A . Since the exact form of this term is 
unknown, we write the exact wavefunctions in the simpler form 

^ix) = e-^fix), (4.3) 

where the functions f{x) must satisfy the condition /(O) = and remain to be 
determined by the iteration method. Substituting (4.3) into (4.2), we obtain 

''^^-2#+fl-E+4 + ^V. (4.4) 



dx"^ dx \ x°^ 



X 
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With Ao(x) — 2x and so{x) = 1 — E + ^ + '^^'J^^'^ , we may compute A„(a;) and 
Sn{x) using (2.14) and the eigenvalues may be calculated by means of the condition 
(2.13). For each iteration, the expression 6 = SnXn-i — Sn-iXn will depend on 
two variables E and x . The eigenvalues E computed by means of 5 = should, 
however, be independent of the choice of x . Actually, this will be the case for most 
iteration sequences. The choice of x can be critical to the speed of the convergence 
to the eigenvalues, as well as for the stability of the process. Although, we don't 
have at the moment a specific method to determine the best initial value of x , we 
may suggest the following approaches. For the spiked harmonic oscillator potential 
V{x) = a;^ + ^ , one choice of x = xq could be the value of xq that minimize 
the potential V , namely xq = (^)^/('^+^) . Another possible choice comes from 
noticing that the ground-state energy of the harmonic oscillator or that of the 
Gol'dman and Krivchenkov Potential can be obtained by setting sq ~ . We may 
therefore start our iteration with xq that is obtained from sg = . For example, if 
a = 4 , then 

A , 7(7 + 1) 



so = l-E+ — + ' = 

X^ X'^ 



implies 



xo = \Ip+\Ip'' + ^—^ (4.5) 



where p = 2(2- 1) ■ "^^^ results of our iteration method for the cases of a = 1.9 and 
ct = 2.1 is reported in Table 1 where we compute the eigenvalues E^^ by means of 
12 iterations only. It is should be clear that these results could be further improved 
by increasing the number of iterations. For the case of a = 4 we report our results 
in Table 2 wherein we used xq to start the iteration procedure, as given by (4.5). 
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Table (I) A comparison between the exact eigenvalues Eni for dimension 
= 2 — 10 of the Hamiltonian (4.2) computed by direct numerical integration of 
Schrodinger's equation and the eigenvalues computed by means of the present 
work. 





a = 1.9 


a = 2.1 


N 


Eoo 


Eqo 


E21 


E21 


2 


8.485 38 


8.485 45 


16.543 63 


16.543 76 


3 


8.564 36 


8.564 42 


16.904 44 


16.904 42 


4 


8.795 44 


8.795 47 


17.381 71 


17.381 45 


5 


9.163 09 


9.163 09 


17.955 44 


17.955 22 


6 


9.646 70 


9.646 68 


18.607 07 


18.607 00 


7 


10.225 04 


10.225 03 


19.320 69 


19.320 73 


8 


10.879 07 


10.879 07 


20.083 41 


20.083 46 


9 


11.592 98 


11.592 98 


20.885 02 


20.885 03 


10 


12.354 18 


12.354 18 


21.717 61 


21.717 59 



Table (II) A comparison between the exact eigenvalues Eq for 1-dimension of 
the Hamiltonian (4.2) with a = 4 computed by direct numerical integration of 
Schrodinger's equation and the eigenvalues Eq computed by means of the present 
work. 



A 


7 


EE 


^0 


0.001 


3 


9.00011427833 


9.00011427912 




4 


11.00006349067 


11.00006349074 




5 


13.00004040373 


13.00004040364 


0.01 


3 


9.00114219619 


9.00114219940 




4 


11.00063478892 


11.00063478889 




5 


13.00040400063 


13.00040400060 


0.1 


3 


9.01136393266 


9.01136402618 




4 


11.00633609974 


11.00633609923 




5 


11.00403643257 


13.00403643252 


1 


3 


9.108660360401 


9.10865860752 




4 


11.06224182608 


11.06224171938 




5 


13.04001518318 


13.04001518306 
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4.2 The Quartic Anharmonic Oscillator Potentials 

We investigate the Schrodinger equation Hip = E{A)'il; for the quartic anharmonic 
oscillators [15]-[19], where 

H = -f^+x^ + Ax'' (4.6) 

In order to obtain the energy levels using the iteration method, we write the exact 
wavefunctions in the form 

ip{x) = e 2 ji^x) 

Consequently, after substituting in the Schrodinger equation, we obtain 

or \q{x) = 2x and so{x) = 1 — E + Ax'^ . We start the iteration in this case with 
xq = 0, the value of x at which the potential takes its minimum value. We report 
our computational results in Table 3. 



Table (III) A comparison between the eigenvalues , n = 0,1,..., 5, for 
the quartic anharmonic oscillator with A = OA computed by direct numerical 
integration of Schrodinger's equation [19] and the eigenvalues E^ computed by 
means of the present work. 



n 


EP 


E 





1.065286 


1.065286 


1 


3.306871 


3.306872 


2 


5.747960 


5.747959 


3 


8.352642 


8.352678 


4 


11.09835 


11.09860 


5 


13.96695 


13.96993 
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5. Conclusion 

One can find a Taylor polynomial approximation about xq for an initial value 
problem by differentiating the equation itself and back substituting to obtain suc- 
cessive values of y^''\xo). This method is perhaps as old as the very notion of a 
differential equation. In this paper we develop a functional iteration method related 
to this general idea and specifically taylored for an important base-class of linear 
equations of the form L{y) = y" — \y' — sy = 0. The iteration is assumed either to 
terminate by the condition Sn/Xn = Sn-i/\n-i = or this condition is imposed, 
as an approximation. After looking at some well-known problems which are exactly 
of this sort, our principle application is to Schrodinger eigen equations. These latter 
problems are converted to the base-type by first factoring their solutions in the 
form '(/^(x) — f{x)y{x), where f{x) is the large- a; asymptotic form, and y sat- 
isfies L{y) = 0. Some aspects of this approach, such as the iteration termination 
condition, the construction of asymptotic forms, and the choice of xq, still await 
more careful mathematical analysis. However, even in its present rudimentary state, 
the method offers an interesting approach to some important problems. This is es- 
pecially so, as it turns out, for problems such as the spiked harmonic oscillator that 
are known to present some profound analytical and numerical difficulties. 
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